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We systematically examine the properties of null geodesics around an electrically charged, asymptotically flat 
black hole in Eddington-inspired Bom-Infeld gravity, varying the electric charge of black hole and the coupling 
constant in the theory. We hnd that the radius of the unstable circular orbit for massless particle decreases 
with the coupling constant, if the value of the electrical charge is fixed. Additionally, we consider the strong 
gravitational lensing around such a black hole. We show that the deflection angle, the position angle of the 
relativistic images, and the magnification due to the light bending in strong gravitational field are quite sensitive 
to the parameters determining the black hole solution. Thus, through the accurate observations associated with 
the strong gravitational lensing, it might be possible to reveal the gravitational theory in a strong field regime. 

PACS numbers: 04.50.Kd, 04.40.Nr, 04.70.-s 


I. INTRODUCTION 

Gravity is a fundamental interaction, and there are many observations and experiments to test the gravitational theories, 
some of which seem to rule out theories alternative to the general relativity. These tests did not show any defects in general 
relativity, proposed by Einstein, while they have been implemented basically in a weak field regime iQl]. On the other hand, 
unlike the tests in a weak field regime, the observations in a strong field regime are quite restricted, which are still not enough 
to verify the gravitational theory in such a strong field regime. That is, the gravitational theory might be modified in the strong 
field regime, which can lead to the phenomena depending on the detail of the gravitational theories. If so, one will be able to 
probe the gravitational theory via the observations associated with compact objects. In practice, there are some suggestions for 
distinguishing the gravitational theories by using astronomical observations (e.g., M). 

So far, many modified theories of gravity have been proposed, where some of them are formulated in order to solve open 
problems in general relativity. Eddington-inspired Born-Infeld theory of gravity (EiBI) is also one of them 101, with which one 
can avoid not only the big bang singularity OSU^, but also the singularity appearing due to the gravitational collapse of dust 
im. EiBI is based on the Eddington gravity lUlI and on the nonlinear electrodynamics of Born and Infeld ifl^ . In EiBI, the 
connection is considered as an independent field as well as the metric, following a Palatini formalism. One of the interesting 
features of EiBI is that the theory in vacuum completely agrees with general relativity. One can observe the deviation of EiBI 
from general relativity with matter, where such deviation becomes significant especially inside neutron stars. That is, one could 
distinguish EiBI from general relativity via the observations of neutron stars itself and/or the phenomena associated with neutron 
stars DSIIMll]- 

The discussions about the neutron stars in EiBI are lively since EiBI was proposed, while the examinations with respect to the 
black holes in EiBI are limited. This might be because the simplest black hole solution in EiBI, i.e., the spherically symmetric 
spacetime in vacuum, agrees with the Schwarzschild black hole in general relativity. Nonetheless, EiBI with non-zero energy- 
momentum tensor can deviate from general relativity even if there is no matter. In other words, en electrically charged black 
hole in EiBI is expected to be different from that in general relativity. Such a black hole solution is discussed in Refs. |Ia I20| - I^ . 

With respect to the phenomena around black holes, the light bending due to strong gravitational field is one of the important 
properties from the observational point of view. In particular, when the light passes in the vicinity of the black hole, the deflection 
angle can become more than 27r, where the light goes around the black hole before approaching the observer. These phenomena 
are known as the strong gravitational lensing. In general relativity, the properties of such phenomena are examined well around 
various black hole solutions ||23| - I^ . In a similar vein, recently the strong gravitational lensing around the electrically charged 
black hole in EiBI is partially examined 1^ . However, such examinations might be insufficient, where only the case of 
the positive coupling constant in EiBI is considered ll^ or the angular position and magnification of the relativistic images due 
to the gravitational lensing are not determined ll^ . Thus, in this paper, we systematically examine the properties of strong 
gravitational lensing around the electrically charged black hole in EiBI, varying the both positive and negative coupling constant 
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in EiBI in addition to the value of the electrically charge. For this purpose, we derive the analytic formulae describing the strong 
gravitational lensing, basically according to Ref. j^ . Then, we will show that the various properties in EiBI are significantly 
sensitive to not only the coupling constant in EiBI but also the electrically charge of black hole. In this paper, we adopt geometric 
units, c = G = 1, where c and G denote the speed of light and the gravitational constant, respectively, and the metric signature 

is (—, +) +, +)■ 


II. ELECTRICALLY CHARGED BLACK HOLE IN EIBI 


EiBI proposed by Banados and Ferreira ||3l is described with the action 

-S'= 4'm], (1) 

where is the symmetric part of the Ricci tensor constructed by the connection E^^, denotes the absolute 

value of the determinant of the matrix {g^^, + kR^^), 9 denotes the determinant of the metric g^,y, Sm denotes the action 
of matter fields, and T'm collectively denotes any matter field. As seen in Eq. (ID, EiBI has two parameters, i.e., k and A. 
K is the Eddington parameter with the dimension of length squared, while A is a dimensionless constant associated with the 
cosmological constant A as A = (A — 1)/k. Additionally, since the action of EiBI in the limit of k = 0 or for = 0 reduces 
to the Einstein-Hilbert action I1 a I^| 2^ . EiBI in the limit of k = 0 and/or Sm = 0 becomes equivalent to general relativity 
with the cosmological constant. The Eddington parameter k is constrained via the solar observations, big bang nucleosynthesis, 
and existence of neutron stars 0 HH-dU] . In addition to these constraints, the possibilities to observationally constrain k are 
suggested with the simultaneous measurements of the stellar radius of the O.SMq neutron star and the neutron skin thickness of 
208pb lllTtl, and with the frequencies of the neutron star oscillations ifl^ . 

Since in EiBI the connection E^^^ is considered as a field independent of the metric the field equations can be obtained 
from the variations of the action Q with respect to the connection and metric J^l; 


Qfiiy — 9^1^ T ^Rfii>: (2) 

, (3) 

where q^i, is an auxiliary metric associated with the connection as E^^ = q^^' {qaa,p + qcTp,a — qap,a) /2, q is the determinant 
of q^^, and denotes the standard energy-momentum tensor. We remark that q^’' is just matrix inverse of q^^, i.e., q^’' ^ 
gf^°‘g'^Pq^p and q^°‘qua = 5^- Meanwhile, raising and lowering indices in T^i, should be done with the physical metric 
The metric describing the spherically symmetric objects is expressed as 

q^^dx^dx'' = -ip'^fdt^ + f-^dr^ + r'^dn^, (4) 

qfj,i,dx^dx'' = —G^Fdt^ + F~^dr^ + Fl^dQ^, (5) 


where tp, /, G, F, and FI are functions of r, and dfl^ = d 6 ^ + sin^ Odcp^. To consider an electrically charged black hole, we 
adopt the energy-momentum tensor given by = [ 9 °"^Ffj^aF^g — 9 ^vFapF°‘^/ /Ett. Then, one can derive the electrically 
charged black hole solution in EiBI as 


/ 


r^Xr'^ F kQ'^ 
Ar^ — kQ"^ 


J \ Ar"^ -I- kQ"^ 


iP = 


'/Xr'^ 


yjXr^ F 




, 0,0 , 


( 6 ) 

(7) 

( 8 ) 


where M and Q denote the mass and electric charge of the black hole, while is the electric field outside the black hole 
13,113,mil- With these metric functions in g^,j, the metric functions in q^i, are determined by 



G = ip 
F[ = r 



(9) 

( 10 ) 

( 11 ) 
















3 


We remark that the black hole solution in the limit of Q = 0 reduces to the Schwarzschild-(anti-) de Sitter spacetime in general 
relativity, i.e.,/(r) = 1 — 2M/r — Ar^/3 and'0(r) = 1, while in the limit of/c = 0 reduces to the Reissner-Nordstr6m-(anti-)de 
Sitter spacetime in general relativity, i.e., /(r) = 1 — 2M/r + jr^ — Ar^/3 and ^(r) = 1. Hereafter, for simplicity, we 
consider the asymptotically flat black hole solution, i.e., A = 1 (A = 0). In this case, there is a maximum value of Q/M 
depending on k/M"^ below which the black hole solution can exist. Such maximum values of Q/M are shown as a function of 
k/M^ in Fig. 5 in Ref. ll^ . 


III. NULL GEODESIC 


We consider the geodesic equation around the electrically charged black hole in EiBI. Since the Lagrangian given by £ = 
g^v{dx^/dT)[dx''/dr)/2 is conserved along the geodesic, i.e., dC/dr = 0, one can put C = —1/2 for a massive particle using 
the scaling degree of freedom of r while C — 0 for a massless particle independently of the scaling degree of freedom of 
r. The motion of particle is subject to the Euler-Lagrange equation, 

_^_d_ f dC 

dxi^ dr \ dx>^ 

On the other hand, adopting the metric ansatz (IDl, we obtain 

2£ = + f-^f^ + r^{e^ + sin^ (13) 

Due to the nature of the static spherically symmetric spacetime, dC/dt = dC/d(/> = 0 in Eq. (fT^ . which leads to that 
the derivatives of t- and (/-components of the four velocity with respect to r become zero, i.e., Pt = P 4 > = 0, where 
= dC/dx^ = . Thus, one can get 



e - 

t = ^ and 0 : 

fV 




(14) 


where e and £ are constants corresponding to the energy and angular momentum, respectively, of the massless particle. In the 
case of a massless particle (£ = 0), rescaling the affine parameter as t —r/e and introducing the impact parameter b = £/e, 
the equation for the radial motion becomes 


^2 



+ 1 / = 1 , 


V{r) 


J.2 ’ 


(15) 


where to derive Eq. ( fTST i we assume that the motion is confined on the equatorial plane (9 = 7r/2). We remark that the orbit of 
massless particle depends only on the value of b independently of the individual values of e and £. 

The unstable circular orbit (UCO) or the photosphere, whose radius is denoted by Ruco, is obtained by solving dV/dr = 0 
for d?V/dr‘^ < 0. Then, the massless particle moves on the UCO around the black hole with infinite time, if U(i?uco) = 1- 
Otherwise, the particle infalling from the infinity either plunges into the black hole if U(i?uco) < 1 or is scattered off by 
the potential barrier at some radius larger than i?uco if y(,R\]Co) > 1- The condition U(i?uco) ^ 1 becomes equivalent 
to 6^ ^ b/., where he is a critical impact parameter depending merely on the black hole parameters. We remark that, for the 
Reissner-Nordstrom case (k = 0), i?uco and be can be obtained analytically as 


3M 


bi = 



(16) 

f^UCO 

f-r)(f?uco 

(17) 


where r± = M ± ^— Q/^. 

As an example, in Eig. [U we show the effective potential V{r) in EiBI with Q/M = 0.5 and k/M^ = 6, where the solid, 
dotted, and broken lines correspond to the cases for b = 2M, be, and 7M, respectively. In the same figure, the position of i?uco 
in EiBI is denoted by with the vertical solid line, where the position of i?uco in general relativity, i.e., for Q/M = 0.5 and 
K = 0, also denoted by the vertical dot-dashed line for reference. Erom this figure, one can see that i?uco in EiBI can deviate 
from that in general relativity. In fact, i?uco = 2.67M in EiBI with Q/M = 0.5 and k/M'^ = 6, while i?uco = 2.82M 
in general relativity, i.e., the deviation is 5.3%. The value of be in EiBI also can deviate from that in general relativity, i.e., 
be = 4.83M in EiBI with Q/M = 0.5 and k/M'^ = 6 while he = 4.97M in general relativity. Additionally, the particle with 
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FIG. 1: Effective potential V{r) in EiBI with Q/M = 0.5 and k/M^ = 6. The solid, dotted, and broken lines correspond to the cases for 
b = 2M, be, and 7M, respectively. The position of Ruco in EiBI is denoted by the vertical solid line, while that in general relativity is 
denoted by the vertical dot-dashed line for reference. The massless particle with b — 7M is scattered by the black hole at r = ro, which is 
denoted by the filled circle. 



QIM Klhf 


FIG. 2: Radius of the unstable circular orbit i?uco around the electrically charged black hole in EIBI as a function of the electrical charge 
Q/M with the fixed value of the coupling constant ac/M^ (left panel) and as a function of ac/M^ with fixed value of Q/M (right panel). In 
the figure, the thick-solid line corresponds to the result in general relativity, i.e., for the Reissner-Nordstrbm spacetime in the left panel and for 
the Schwarzschild spacetime in the right panel. 


b > be is scattered by the black hole, as mentioned the above, i.e., for instance, the particle with b = 7M is scattered at r = ro, 
which is also denoted by the filled circle in Fig. [T] 

In Fig. |2] we show the dependence of i?uco in EiBI on the electrical charge Q/M with fixed value of the coupling constant 
k/M'^ in the left panel and on k/M"^ with fixed value of Q/M in the right panel, where the labels in the figure denote the 
fixed value of k/M"^ in the left panel and the fixed value of Q/M in the right panel. For comparison, we also show i?uco in 
general relativity with the thick-solid line, i.e., for the Reissner-Nordstrom spacetime in the left panel and for the Schwarzschild 
spacetime in the right panel. From this figure, one can observe that i?uco depends strongly not only on the value of Q/M but 
also on the value of k/M^. In particular, i?uco with a fixed value of Q is a monotonically decreasing function of the coupling 
constant k. Note that the unstable circular orbit exists, just like the event horizon even for Q/M > 1, provided ac < 0. On 
the other hand, the critical impact parameter be with which the massless particle moves on the UCO, also depends on the values 
of Q/M and k/M'^. In Fig. [3 we show the dependence of 5c on Q/M in the left panel and on k/M"^ in the right panel, as 
same as in Fig. |2] From this figure, we find that the dependence of be on Q/M and k/M'^ are qualitatively the same as that of 
f?uco- We also find that be in EiBI cannot be over be for the Schwarzschild spacetime, which is be = S'/SM, in the range of 
the parameters we adopt in this paper. 
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FIG. 3: Critical impact parameter be in EiBI as a function of the electrical charge QjM with the fixed value of the coupling constant k/M^ 
(left panel) and as a function of k/M^ with fixed value of QjM (right panel). As same as in Fig|^ the thick-solid line corresponds to the 
result in general relativity. 


IV. DEFLECTION ANGLE 


We especially consider the scattering orbit of massless particle around the electrically charged black hole in EiBI, i.e, for 
b > be. Assuming the particle motion on the equatorial plane and using Eqs. (fT4l) and dTSl ). one can derive the equation, 


d(j) 

dr 


d(j)/d7 


bip 


(18) 


dr / dr — fip'^b^ 

As shown in Eig. IH we consider that the particle turns around the black hole at r = ro and (p = 0. We remark that the position 
of ro in the effective potential (Eig. [T]i corresponds to this turn-around point. Since, by definition, dr/dp should be zero at the 
turn-around point, the impact parameter for b>bc can be expressed as 


V = 


(19) 


/(ro)V'^(ro)' 

In Eig. |5] the impact parameter determined by Eq. (fT9l l is shown as a function of the position of the turn-around point rg, 
where the left and right panels correspond to the results for QjM = 0.5 and 1.0, respectively. In the figure, the labels 
denote the values of k/M^, where the case in general relativity (n = 0) is also plotted by the thick-solid line for refer¬ 
ence. The leftmost end-point of each curve corresponds to the critical impact parameter be for each case, which is given by 
be = f?uco/[/('^uco)V^^(^uco)]- Erom this figure, we find that the impact parameter in EiBI for the massless particle scat¬ 
tering far from the black hole, such as tq > 5M, is almost independent from the coupling constant, while that depends strongly 
on the coupling constant and the electrical charge for the scattering orbit close to the black hole. 

Additionally, integrating Eq. (fT^ . we can get the equation. 


p{r) - p{ro) = f 

J rn 


btjj 


r-\/r^ — fip'^b'^ 

Therefore, as shown in Eig. |4l the deflection angle A(p{ro) is determined as 


dr. 


( 20 ) 


A(p = 2p(oO) — TT 


= 2b 




rp'r^ - 


dr — TT, 


( 21 ) 


where we adopt pjrp) = 0. However, one can see that the integrand in Eq. (1211 1 diverges at r = rp due to Eq. (fT9l l. That is, as 
shown in 11^ 1^ , one should remove the pole at r = tq for determining the deflection angle A(p. 

Now, introducing a new variable z defined as z = 1 — rp/r, i.e., r(z) = ro/(l — z), the angle p{oo) can be express as 


where is given by 


I(ro) = p{oo) = [ p{z)T{z,ro)dz, 
Jo 

br{z) 


( 22 ) 


^{z,ro) = 


ro\/r'^{z) - f{z)p‘^{z)b'^ 


(23) 
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FIG. 4: Image of scattering orbit of massless particle, where Aip is deflection angle. The right panel is magnified view in the vicinity of the 
scattering position, where ro is the curtate distance between the particle orbit and the black hole. 



FIG. 5: For the scattering massless particle, impact parameter in EiBI is shown as a function of the position of the tum-around point, where the 
left and right panels correspond to the results for Q/M = 0.5 and 1.0, respectively. The labels in the figure denote the value of the coupling 
constant For reference, the case in general relativity (k = 0) is also shown with the thick-solid line. 


As mentioned the above, the function of ro) diverges at z = 0 (r = ro), while t/’(z) is regular for any values of z. To see 
the behavior of divergence of the function at z = 0, we approximate it near z = 0 as 


^{z,ro) - , = To{z,ro), 

yjaz + /3z2 


(24) 


where a and /3 are functions of rp given as 

a(?'o) 


P{ro) 


V dz \ ^=0 

2/ot/’o - roiMoY^ 
-/ot/’o^+ro(/oV'o7-y(/oV'o")"- 


(25) 


(26) 


Here, the variables with the subscript 0 denotes the corresponding values at r = tq or z = 0, and the prime denotes the 
derivative with respect to vq. Thus, one can find that the value of T(ro) is finite when a is nonzero. Meanwhile, as mentioned in 
the previous section, dV/dr = 0 is realized only at r = R\jco^ which leads to the statement that a = 0 only at tq = i?uco ■ So, 
the deflection angle is finite value for tq > i?uco and diverges at ro = i?uco ■ 

In order to determine the deflection angle by removing the singular point, we separate T(ro) into two parts as 


d^(ro) = 'd^D{ro) + 2 :fl(ro), 


( 27 ) 
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FIG. 6: Deflection angle is shown as a function of the position of the turn-around point. The left and right panels correspond to the results for 
QjM = 0.5 and 1.0, respectively. The labels in the figure denote the value of the coupling constant ac/M^ in EiBI, where the case in general 
relativity (k = 0) is also shown with the thick-solid line. 


where 


^D(ro)= / 'il^QFQ{z,rQ)dz, 

JO 

(28) 

1-R{ro)= / g{z,ro)dz, 

(29) 

ro) = i/{z)R{z, ro) - i/oTo{z, tq). 

(30) 

The first term in Eq. dZTli can be analytically integrated which becomes 


T r \ 2V’o, y/fd + y/a + fd 
ddD{ro)= ^log /- 

V? Va 

(31) 


On the other hand, for tq > Rvco^ the function Q{z, tq) can be expanded around z = 0 as 


G(z,ro) = [ipiz) - 'ij}o]To{z,ro) 

d^{z) 


dz 


«(fo) 


Oiz), 


(32) 


which leads to Q{0, tq) = 0, i.e., I/{(ro) becomes a finite value. Finally, one can calculate the deflection angle for tq > i?uco 
via the relation. 


A(p{ro) = 2lD{ro) + 2XR{ro) - tt. (33) 

In Fig. |6] we show the deflection angle as a function of ro, where the left and right panels correspond to the results with 
Q/M — 0.5 and 1.0, respectively. The labels in the figure denote the adopted values of the coupling constant k/M'^ in EiBI. 
In any case, the deflection angle increases as ro decreases, and becomes equal to 27r with a specific value of ro, where the 
massless particle goes around the lens object before reaching the observer. As ro additionally decreases, the deflection angle 
increases more and the number that the massless particle goes around the lens object also increases. When ro eventually reaches 
the position of the unstable circular orbit, the deflection angle diverges, where the massless particle goes around the lens object 
for all eternity. The positions of the relativistic images where the deflection angle becomes equal to a multiple of 27r depend 
strongly on the gravitational theory, as seen in Figl6] In fact, the positions of the relativistic images for the negative (positive) 
coupling constant in EiBI become larger (smaller) than that expected in general relativity, where the difference between the 
positions expected in EiBI and in general relativity seems to increase with Q/M. On the other hand, since the deflection angle 
for ro ^ 5M is almost independent from the gravitational theory, it could be quite difficult to distinguish the gravitational theory 
via the observation of the deflection angle for rg > 5M even if detected. 
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V. STRONG GRAVITATIONAL LENSING 

Here, we focus on the deflection angle for rg ~ i?uco^ where a(ro) is expanded as 

a = ai(i?uco)(’’o - R\JCo) + C>((^o - -Ruco)^), 
ai(?'o) = ifoipo)' - roifoipo)". 

Substituting this expansion into Eq. (ISTT l. one can get 


loiro) = -ai log 


^■o 


,^uco 

where oi and 02 are the constants determined at rg = Rtjco as 

at = '0o/V/3(i?uco), 
4/3(i?uco) 


— 1 I +0,2 + 0(rg — i?UCo), 


02 = Oi log ■ 


^ucoai(i?uco) 

The second term in Eq. dZTl i can be expanded in the vicinity of rg = i?uco as 

nl Qng 

—ATg - rtuCOf / 

nl I, 

n—0 


00 

Rniro) = “ RvcoT / 

Jo 


drJl 


0 ^’0 


dz 


i'o=Rvco 


f Q{z,R\jco)dz + 0 {rg — Ruco)- 
Jo 


Taking into account that Q;(i?uco) = 0, G{z, i?uco) can be expressed in the vicinity of z = 0 as 


G{z, Rtjco) = [ij{z) -'J’o]J^o{z, Rtjco) 


d^l^iz) 


dz 


1 


!=o \//3(i?uco) 


0{z). 


That is, Xji{rg) is still a hnite value. Thus, one can calculate the deflection angle with 

ro 


Aip{ro) = -2oi log 


— 1 ) + 02 + O(ro — i?uco), 


-Ruco 

02 = 25,2 + 2Xji[RTjco) — 

Eurthermore, from Eq. (fT^ . one can get the expansion of h{ro) around ro — Rtjco as 




(ro - -Ruco)^ ■ 


ro=Rgco 


With this expansion, the deflection angle [Eq. (HTt l is rewritten as 

A(p{b) = -oi log + 03 + 0{{b - 

03 = O 2 + Oi log 


i?2 U2 


wco 


d^b 


2bc drl 


ro=Rgco , 


(34) 

(35) 


(36) 

(37) 

(38) 


(39) 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 


Now, we consider the lens geometry shown in Fig. |7] We assume that the observer and source are located in the flat spacetime, 
while the region around the black hole is strongly curved spacetime, which leads to the light bending. That is, the light ray passes 
from the point S up to the point O via the point C in Fig. |7] Then, the lens equation can be expressed as 


tanw = tan©-[tan© + tan {Aip — ©)] 

Dos 


(46) 


where uj corresponds to the angular separation between the lens and source, and © corresponds to that between the lens and 
image, i.e., ui = ZEOS and © = ZLOI in Fig. |7] One can also see that b = 79ol sin ©. 
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FIG. 7: Lens diagram, where S, L, O, and I are the positions of the source, lens object, observer, and image, respectively. 



FIG. 8: The angle of the 1st relativistic image, 0?, in EiBI as a function of the electrical charge Q/M with the fixed value of the coupling 
constant KjM^ in the left panel and as a function of k/M^ with the fixed value of QjM in the right panel. The labels in the figure denote the 
fixed values of Q/M or k/M/^. For reference, the solid-thick line in each panel corresponds to the result in general relativity. 


In the strong deflection limit, i.e., w <C 1, 0 -C 1, and (A(/3„ — 0) <C 1, the lens equation reduces to 

= 0 - (47) 

Mos 

where A<p„ is the deflection angle removed all the loops of photon around the lens object ll^ . That is, A<p„ = A(p — 2mr with 
an integer n for 0 < A(p„ <C 1. In this limit, since b ~ Z7ol 0, the deflection angle [Eq. (I44l i1 reduces to 


A<^(0) 


-oi log 


/ Z7ol0 

V b. 



+ as- 


Using Eq. (|48]) . one can get the angle 0°, with which A(p is equal to 2mT, i.e., 

be 


= 


e„ = exp 


'as — 2mr 


ai 


(48) 


(49) 

(50) 


In Eig. m we show the angle of the 1st relativistic image, 0°, determined from Eq. ( |49] | as a function of the electrical charge 
Q/M with the fixed value of the coupling constant in EiBI k/M"^ in the left panel and as a function of the value of k/M'^ with 
the fixed value of Q/M in the right panel, where we assume that Dql = 8.5 kpc and M = 4.4 x 10 ®Mq to match the previous 
calculations in 12^ . Erom this figure, we And that the angle of relativistic image for the fixed value of nonzero Q/M decreases 
as the coupling constant increases. In particular, the deviation from the result in general relativity significantly increases with 
the value of |k/M^ | and Q/M. This is consistent with the result in Eig. |6l where the position of the turn-around point decreases 
as the coupling constant increases. In fact, compared with the results in general relativity, the value of 0® with Q/M = 0.5 in 
EiBI becomes 3.4% larger for = —10 and 4.8% smaller for = +10, while 0® with Q/M — 1.0 in EiBI becomes 

22.5% larger for = -10. 
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Setting 0 = 0° + A0„ (| A0„/0° | <C 1) and expanding A(p{Q) around 0°, one can derive Aipn as a function of A0„ 


_ aiHoL.^ 

— 7 ^©71 


(51) 


Then, the lens equation becomes 


iO — + ( 1 + 


ai-PpL T>ls 
^OS 

^ qO,^i£olD^^ 

— 7 ^ 7-j 

T/qS 


A 0 , 


(52) 


where we assume Z?ol ^ to obtain the final expression. Then, the position of the n-th relativistic images, 0 = 0„, are given 
by 


,0 I bcCnDoS 

'n — 'Jji 


e. = 0 ^,+ (^-e°). 


aiD-LsDoi, 

Additionally, the magnification of the images defined by the inverse of the Jacobian estimated at the images, /x„, is given by 


(53) 


Mn — 


w duj 

©M 


e=eo 


bpenPos 
oiZJol^'ls ’ 


(54) 


where we simply consider the image at 0 = 0 ° and we adopt the relation of ZJql ^ be again ll^ . 

Finally, according to Ref. we consider the simplest situation for observation. That is, we assume that only the outermost 
relativistic image 0^ can be resolved from the others and all the remaining images are concentrated at 0j^. Then, the observable 
quantities are 


s = 0? - 0[J^ = 0^ exp(a 3 /ai - 27r/ai), 


-1 


11= X! 

\ ra =2 


~ exp(27r/ai), 


(55) 

(56) 


where we adopt the relations of exp(27r/ai) ^ 1 and exp(a 3 /ai) ~ 0{1) to derive the right hand side of Eq. (l56l l. We remark 
that the value of TZ is dependent on the properties of black hole spacetime, and independent of Dos, ^ls, ^OL, and oj. We also 
remake that the values of oi and 03 can be determined from the observable quantities s and TZ as an inverse problem. 

For example, we consider a specific case with the supermassive black hole located at the center of our Galaxy, whose mass 
is estimated to be M = 4.4 x 10®M© 1^ . while the distance from the solar system, ZJol, is around 8.5 kpc 1^ . In Fig. 
|9] we show the angular separation between 0® and 0^ as a function of the electric charge Q/M with the fixed value of the 
coupling constant k/M'^ in the left panel and as a function of the value of k/M'^ with the fixed value of Q/M in the right 
panel. From the right panel of this figure, one can observe that the angular separation s with the fixed value of nonzero Q/M 
increases with the coupling constant. In most cases, the separation angle in EiBI becomes larger than that for the Schwarzschild 
black hole, but s can be smaller for k/M"^ < —5 and Q/M < 1.5. Compared with the results in general relativity, i.e., for the 
Reissner-Nordstrom spacetime, the angular separation s with Q/M = 0.5 in EiBI becomes 27.0% smaller for /t/M^ = —10 and 
92.6% larger for n/M'^ = +10, while s with Q/M = 1.0 in EiBI becomes 82.5% smaller for n/M'^ = —10, which could be a 
significant deviation from the expectation in general relativity. Moreover, in Fig. (TO] the magnitude of the relative magnification 
TZm, defined as 


= 2.51oglQ7^, (57) 

is plotted as a function of the electric charge Q/M with the fixed value of the coupling constant k/M"^ in the left panel and as a 
function of the value of k/M'^ with the fixed value of Q/M in the right panel. We find that, with the fixed value of nonzero Q/M, 
the relative magnification decreases as the coupling constant increases in any case. In particular, in the case with k/M^ < —6 
and Q/M < 1.3, the relative magnification can be larger than that for the Schwarzschild black hole, which is the maximum value 
predicted in general relativity. That is, if the relative magnification larger than that expected for the Schwarzschild black hole 
would be observed, it could be the observation detecting an imprint of an alternative gravitational theory, not general relativity. 
Moreover, compared with the results in general relativity, the magnitude of the relative magnification with Q/M = 0.5 in EiBI 
becomes 5.5% larger for k/M^ = —10 and 12.7% smaller for k/M^ = +10, while TZm with Q/M = 1.0 in EiBI becomes 
42.5% larger for k/M^ = —10. 
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FIG. 9: Angular separation s in EiBI as a function of the electric charge Q/M with the fixed value of the coupling constant in the left 

panel and as a function of the value of k/M^ with the fixed value of Q/M in the right panel. 



FIG. 10: Magnitude of the relative magnification TZm as a function of the electric charge Q/M with the fixed value of the coupling constant 
k/M^ in the left panel and as a function of the value of ac/M^ with the fixed value of Q/M in the right panel. 


VI. CONCLUSION 

In this paper, we systematically examine the properties null geodesics around an electrically charged black hole in Eddington- 
inspired Born-Infeld gravity (EiBI), where we especially focus on the asymptotically flat black hole solution. In this case, the 
free parameters determining the back hole solution in EiBI are the coupling constant in EiBI (k/M^) and the electrical charge 
{Q/M), which are normalized by the mass of black hole M. In order to examine the orbital motion of a massless particle, we 
derive the null geodesic. We And that the radius of the unstable circular orbit, i?uco^ depends strongly on both and Q/M, 

where i?uco with the fixed value of nonzero Q/M decreases as k/M^ increases. 

We also examine the strong gravitational lensing around the electrically charged black hole in EiBI. We show that the de¬ 
flection angle due to the light bending in strong gravitational held diverges, when the distance between the black hole and the 
position of the turn-around point of light agrees with i?uco as same as in general relativity. The position angles of the relativistic 
images also strongly depend on both k/M^ and Q/M, where the angle with the fixed value of nonzero Q/M decreases with 
k/M"^. Eurthermore, in the strong deflection limit, we derive the analytic formulae to determine the angular position and magni¬ 
fication due to the winding of light around the black hole in EiBI, according to Ref. ll^ . In the simplest situation where only the 
outermost image can be resolved from the others, we concretely calculate the angular separation and the relative magnification, 
supposing the central black hole in our Galaxy. Then, we And that the angular separation and the relative magnification are quite 
sensitive to k/M^ and Q/M. In particular, the relative magnification in EiBI for k/M^ < —6 and Q/M < 1.3 can be larger 
than the maximum value expected in general relativity. Since these observable quantities are directly affected by the parameters 
determining the black hole solution, it might be possible to make constraints on such parameters via the observation of the strong 
gravitational lensing ImI] . Even if the current astronomical instruments are not enough to test the gravitational theory in a strong 
held regime, these observational constraints could enable us to reveal the gravitational theory in the future. 
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